Fermi surface of SrFe2P2 determined by de Haas-van Alphen effect 
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We report measurements of the Fermi surface (FS) of the ternary iron-phosphide SrFe2P2 using 
the de Haas- van Alphen effect. The calculated FS of this compound is very similar to SrFe2As2, the 
parent compound of the high temperature superconductors. Our data show that the Fermi surface is 
composed of two electron and two hole sheets in agreement with bandstructure calculations. Several 
of the sheets show strong c-axis warping emphasizing the importance of three-dimensionality in the 
non-magnetic state of the ternary pnictides. We find that the electron and hole pockets have a 
different topology, implying that this material does not satisfy a (tt, tt) nesting condition. 



Many theories of the superconductivity and magnetism 
in the iron-pnictides have their unusual Fermi surface 
topology as a central ingredient [H, [2, H, 01 • Both super- 
conductivity and magnetism can be enhanced by geomet- 
rical nesting of the hole and electron Fermi surface sheets. 
Bandstructure calculations suggest that 1111-arsenides 
(ReFeAsOi-ajFa;, Re e La^Sm) have stronger 'nesting' 
peaks in the non-interacting susceptibility than their 122 
arsenide (XFe2As2, X S Eu, Ba, Sr, Ca) counterparts 
and the former generally have a higher Tc- The 122 ma- 
terials are calculated to have stronger warping of their 
quasi-two-dimensional sheets and experimentally this is 
reflected in a more isotropic upper critical field Hc2- 

One argument against a predominant role of Fermi sur- 
face nesting is that the phosphide analogues of the 122 ar- 
senide materials (XFe2P2) do not show superconducting 
or magnetic order despite their calculated Fermi surfaces 
being very similar to the corresponding arsenides. How- 
ever, partial substitution of As by P in EuFe2(Asi_a;Px)2 
and BaFe2(Asi_xPK)2 suppresses the antiferro- 
magnetic order and results in superconductivity with Tc 
up to 26 K and 30 K respectively. The evolution from an- 
tiferromagetism to superconductivity and then to para- 
magnetism in this series provides a strong test of our 
understanding of the physics of these materials. A par- 
ticular question is whether this is driven by changes in 
the Fermi surface topology. Although bandstructure cal- 
culations can provide a guide, accurate experimental de- 
terminations are essential because inaccuracies in the cal- 
culated band energies (often of order 0.1 eV 0]) can lead 
to significant changes in the topology. 

Quantum oscillation (QO) studies prove a direct way 
to probe in detail the bulk, full three dimensional Fermi 
surface, and can also resolve the strength of the many- 
body interactions at the Fermi level. However, so far QO 
experiments have not been possible on the high Tc super- 
conducting arsenides because of their high Hc2 and dis- 



order induced by doping. Here we present a de Haas- van 
Alphen effect determination of the Fermi surface topol- 
ogy and effective masses of SrFe2P2 which is the end 
member of the SrFe2(Asi_a;Pa;)2 series. Wc find that the 
hole and electron sheets have quite different topology im- 
plying that a geometric nesting condition is not satisfied. 
This perhaps, explains the lack of either superconduc- 
tivity or magnetic order in this compound. The many 
body mass renormalisation are strongly sheet dependent 
and are much larger than expected from simple electron- 
phonon coupling. 

High quality single crystals of SrFe2P2 with residual 
resistivity ratios p(300K)//3(1.8K) greater than 50, were 
grown from a Sn flux. Torque magnetometry was per- 
formed using piezoresistive microcantilevers in high fields 
Q . Band structure calculations were performed using an 
augmented plane wave plus local orbital method as imple- 
mented in the WIEN2K code ^. For these calculation 
we used the crystallographic parameters a = 3.825 A, 
c = 11.612 A and zp — 0.3521 as determined by Mewis 

Fig. [IJa) shows the raw torque signal for a single field 
sweep up to 18 T. de Hass-van Alphen oscillations are 
observed for fields above ^ 4 T. In Fig. [T] (b) we show the 
fast Fourier transform (FFT) spectrum of the signal (in 
inverse field) measured at two angles 9 = 8° and 6 = 49°, 
taken relative to the c-axis of the crystal and rotating 
towards the [100] direction. At all angles the spectrum is 
dominated by a (split) single peak which we label ai^2, 
and in order to make the other smaller amplitude peaks 
visible we show the square-root of the FFT amplitude in 
the figure. In total seven frequencies are observed, which 
we denote as ai, q;2, /32, 7, <5 and e, and four others are 
harmonics of ai^2- At higher angle additional splitting of 
these peaks is observed which likely comes from a small 
^ 1° mosaic distribution of the crystal orientation. Note 
that the FFT peaks below ^ 0.2 kT are noise related. 
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FIG. 1: dHvA oscillation data in SrFe2P2 at two different 
angles, (a) Raw torque signal with smooth polynomial back- 
ground subtracted (b) Fast Fourier transforms. Note that the 
y-axis is the square-root of the FFT amplitude. 



Each of these frequencies are related to a extremal cross 
sectional area of the FS in momentum space Ak via the 
Onsager relation F={h/27re)Ak- 

The calculated Fermi surface for SrFe2P2 consists of 
two concentric electron cylinders at Brillouin zone corner, 
and two hole cylinders centered at the F point (see Fig. 
[5]) which is very similar to that calculated for the non- 
magnetic phase of SrFe2As2 All of the sheets are 
significantly warped. 

Rotating away from the c-axis we can trace the full 



three-dimensional topology of the FS by looking at the 
evolution of the dHvA frequencies as the cyclotron orbits 
traverse different parts of the FS. Figure [2] (left) shows 
the angle dependence of each of the observed frequencies 
with the field being rotated from [001] to [110] on the 
left and [GDI] to [100] on the right. There is a high level 
of correspondence between the observed frequencies and 
the band structure predictions shown in Fig. O (right). 
The inner electron pocket (band 4) corresponds well to 
the measured a frequencies, both in absolute cross sec- 
tional area (dHvA frequency) and angle dependence (c- 
axis warping). Similarly, band 3 matches the P orbits. 
Note the marked differences in angle dependence shown 
by the /3 orbits when rotated towards either the [100] or 
[110] directions. This is well reproduced in the calcula- 
tions and originates from the highly non-circular shape of 



this Fermi surface sheet (see Fig. [2]). Comparing in detail 
the absolute frequencies of the a and (3 with the calcu- 
lations shows them both to be slightly smaller. A shift 
of the band energies up by 59meV and 49meV respec- 
tively for band 3 and 4 brings them into almost perfect 
alignment (see Fig. [3]) . 

The e orbit agrees very well with the maximum fre- 
quency from the large hole sheet (band 2), but is observed 
over only a narrow range of angle (42-58°). The rea- 
son for this is apparent by looking at the slices through 
the FS shown in Fig. [U where it can be seen that for 
this angle there is a large degree of nesting between the 
warped sections of this band. Indeed, the bandstructure 
calculations show a large increase in the curvature fac- 
tor |c?^Afc/9fcy ]~ 2 (and hence signal amplitude) for these 
angles. 

This leaves frequencies 7 and 6 to assign. At first 
glance it appears as though the minima of either band 
2 or band 1 could account for these. However, by plot- 
ting FcosO versus angle (see Fig. [3]) it is apparent that 
7 cannot be a minimal extremal crosssection because it 
decreases as the polar angle increases. Hence, as shown 
in Fig. [3] the minimum of band 2 has the wrong curva- 
ture to explain this orbit. We therefore conclude that 
this orbit must originate from the maximum of band 1 
as this is the only unassigned orbit with the correct lo- 
cal topology. For the frequencies to correspond we need 
to shift the energies of band 1 down by 110 meV. With 
this shift the curvature is in good agreement with data 
(see Fig. [3]). Importantly, as shown in Fig. |4] this shift 
causes the band 1 FS to be a closed three-dimensional 
ellipsoid. With this assignment the remaining 6 orbit 
must come from the minimal (tubular section) of band 2. 
Only a small change in the warping of this sheet would 
be needed to get complete agreement. 

Fig. [3] shows a comparison of all the observed frequen- 
cies with the shift band structure calculations and Fig. [4] 
shows the resulting changes in the Fermi surface topol- 
ogy. The shifts can be thought of as a fine tuning of the 
bandstructure calculations. For example, for band 3 the 
shift corresponds to a decrease in its volume by 0.044 
electrons per units cell. A test of the consistency of the 
shift is to check for charge neutrality. SrFe2P2 is a com- 
pensated metal so the volumes of the electron and hole 
sheets should be exactly equal. With the shifts bands 1- 
4 contain -0.026, -0.284, 0.197, and 0.108 holes/electrons 
respectively, and so the charge balance remains almost 
exact. The imbalance of 0.005 extra holes per unit cell is 
equivalent to < 2% of the total volume of the hole sheets 
showing the accuracy of our Fermi surface determination. 

The effective masses, extracted [1] by fitting the ob- 
served temperature dependent amplitude of the dHvA 
oscillations to the conventional Lifshitz-Kosevich formula 
[1^ , are compared in Table I to the corresponding band- 
structure values. The mass enhancements are strongly 
sheet dependent and range from A = 1 — m* /rub) = 0.3 
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FIG. 2: (color online). Angle dependence of the dHvA frequencies for SrFe2P2. The left panel of the angle plot shows the 
experimentally observed frequencies as the magnetic field is rotated from [110]—* [001] [100]. The right panel shows the 
corresponding predictions of the band structure calculation. At far right we show the Fermi surfaces associated with each band. 
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FIG. 3: (color online). Observed dHvA frequencies multiplied 
by cos 9 along with band structure predictions for shifted en- 
ergy bands. The upper panel shows the data for the 7 orbit 
on an expanded scale 



for the smallest hole orbit to A — 1.1 for the inner electron 
sheet. These enhancements are similar to those observed 
in LaFePO and are much larger than expected from 
electron-phonon coupling alone [13|](Aep — 0.25). The ta- 
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FIG. 4: (color online). Two dimensional sections through the 
FS in the (110) plane before and after shifting the bands to 
agree with our dHvA measurements. The Brillouin zone edges 
are marked. 



ble also shows the orbit specific dHvA mean-free-paths. 
It is interesting to note that, as for LaFePO 7], the elec- 
tron sheets have the longest mean-free-paths and give the 
strongest dHvA signals. It is not clear to us if there is a 
fundamental reason for this. 

The above illustrates that SrFe2P2 has a FS which is 
highly dispersive in the c-axis and the electron and hole 
pockets are far from fulfilling a nesting condition. This 
may explain the observation of an almost isotropic Hc2 in 
the analogous hole and electron doped superconducting 
122 arsenides 



^ ^ In LaFePO, the hole and electron 

FS sheets are much closer in shape and size 0, [3] and 
this proximity to (tt, tt) nesting may well be a factor in 
determining why LaFePO superconducts while SrFe2P2 
does not. 

The differences with SrFe2As2 appear more compli- 
cated to answer because the phosphides and arsenides 
have much in common in their local structure: the As and 
P ions are isoelectronic, they have a similar atomic radius 
(100pm and 115pm respectively), comparable ionization 
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TABLE I: Measured dHvA frequencies, effective masses (m*) 
and mean free paths {£) , along with the values from the band 
structure calculations. The experimental masses were deter- 
mined at 6^ = 9° expect for the e orbit which is at S = 49°. 
The band structure values are all quoted at 6 = 0° {H\\ [001]). 
The final column show the ratio of measured effective mass 
to the band mass at the same angle (at 9 = 49° mi, — 1.98me 

for 2inax). 
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energies and coulomb interactions [17[ and the Pn-Pn dis- 
tance is large compared with their respective molecular 
bond lengths (except for CaFe2P2, which is more analo- 
gous to the collapsed tetragonal phase of CaFe9As9[l8t). 
However, it is important to notice that the 122 arsenides 
have Pn-Fe-Pn (Pn= pnictide) bond angles close to the 
ideal tetrahedral angle of 109.47°. In the phosphides 
however, these bond angles are much larger, and for 
SrFe2P2 is 116.34° [11]. The Fe-Pn distance is also 
smaller by ~ 0.13 A in the arsenides, leading to a dif- 
ferent bonding overlap. In particular, this results in dif- 
ferences in the calculated pnictogen density of states [ioj] 
which Yildirim has argued can dramatically alter the lo- 
cal magnetism of the Fe atoms 21]. 

The Pn-Fe-Pn bond angle departs from an ideal tetra- 
hedral angle in both LaFePO (119.3°) and SrFe2P2 
(116.34°) in a similar manner and this may be why nei- 
ther of these compounds exhibit magnetic order. This 
naturally accounts for the observed suppression of mag- 
netism as SrFe2As2 is doped with P. However, super- 
conductivity arises at intermediate doping, and this im- 
plies that the magnetic order remains important. We 
suggest that local interactions and a proximity to nest- 
ing work together to enhance Tc- The lower Tc of the 
122 compounds as compared to the 1111 compounds is 
then explained by the absence of nesting [J, [3, 13, 
and the phosphides have a much lower Tc than the ar- 
senides because of the absence of local magnetic interac- 
tions 



22,l23|,l24|,l25l 



In conclusion, we have mapped out the full three di- 
mensional Fermi surface of SrFe2P2, the isoelectronic sis- 
ter compound to the antiferromagnet SrFe2As2. The 
Fermi surface is in good overall agreement with the pre- 
diction of our bandstructure calculations with small shifts 
in the band energies. Unlike LaFePO the Fermi sur- 
face of SrFe2P2 is far from fulfilling a geometric nest- 
ing condition, being composed of warped 2D electron 



cylinders, a strongly warped outer hole cylinder and one 
closed hole pocket. The non-magnetically ordered, non- 
superconducting ground state of SrFe2P2 is likely related 
to the combined absence of nesting and local moment 
interactions, compared to the analogous arsenides. 
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